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■ Abstract. We discuss possible definitions for a slip velocity that describes the 

relative motion between large particles and a turbulent flow. This definition is 
necessary because the slip velocity used in the standard drag model fails when particle 

■ size falls within the inertial subrange of ambient turbulence. By decomposing drag 
<4_H . into steady and stochastic components, we obtain both a steady slip velocity and 

c/2 1 a fluctuating slip velocity. From this decomposition, we propose two definitions for 

stochastic slip velocity. These definitions were selected in part due to their simplicity: 



they do not require filtration of the fluid phase velocity field, nor do they require 
^ (— i the construction of conditional averages on particle locations. A key benefit of this 

O f simplicity is that the stochastic slip velocity proposed here can be calculated equally 

well for laboratory, field, and numerical experiments. The stochastic slip velocity allows 
the definition of a Reynolds number that should indicate whether large particles in 
k*" . turbulent flow behave (a) as passive tracers; (b) as a linear filter of the velocity field; 

or (c) as a nonlinear filter to the velocity field. We calculate the value of stochastic 
slip for ellipsoidal and spherical particles (the size of the Taylor microscale) measured 
C \ in laboratory homogeneous isotropic turbulence. The resulting Reynolds number is 

■ significantly higher than 1 for both particle shapes, and velocity statistics show that 

. particle motion is a complex non-linear function of the fluid velocity. We further 

investigate the nonlinear relationship by comparing the probability distribution of 
fluctuating velocities for particle and fluid phases. 

^ ■ 
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1. Introduction 

It is well known that particles suspended in a turbulent flow can significantly modulate 
the turbulent dynamics even at moderate concentrations [B 121 E]. The ability to 
predict turbulence modulation by particles (change in Turbulent Kinetic Energy (TKE) , 
dissipation, etc.), as well as the dynamics of the suspended particles (aggregation- 
dispersion mechanisms, encounter rates, etc.) strongly depends on our ability to predict 
the forces that the individual particles experience in turbulent flows. 

However, quantifying the forces mediating fluid-particle coupling is nontrivial. 
Neither the stochastic nor deterministic description of turbulent suspensions is currently 
at the level required for useful engineering predictions. The limits of current knowledge 
were summarized eloquently by Qureshi et al. [4] : "[...] writing (not to mention solving) 
the particle equation of motion in the most general case remains a challenge, and only 
limited cases are treated at present." 

The forces which explicitly couple fluid and particle phases are the shear and normal 
stresses felt at the particle-fluid boundary. Many prediction schemes attempt to simplify 
the description by replacing the integral of surface forces with a local body force, the 
familiar drag force Fjj. This method works well when particles move through a steady 
flow with a predictable wake structure. However, drag models lose their predictive 
ability when the particle wake is strongly influenced by stochasticity in the fluid phase 
(most notably fluid turbulence). New models are needed, as discussed below. When 
formulating these new models, it is common to use ideas and formulas carried over from 
the simpler flows in which the standard drag model works well. One such idea is slip 
velocity, which is clearly defined in the simple drag model, but becomes problematic in 
turbulent flows. Our goal in this paper is to specify a stochastic slip velocity in a general 
and unambiguous manner. Doing so will advance the community's ability to formulate 
drag models predicting particle-fluid coupling. 

2. Background 

The stokes drag for a sphere in a creeping flow is: 



where d p is particle diameter and U s is the slip velocity U s = Uf — U p , where Uf is 
an appropriate fluid velocity and U p is an appropriate particle velocity. For equation 
([1]), U s is defined as the difference between particle center-of-mass velocity u p and the 
uniform velocity found outside the particle Uoo. When the particle Reynolds number 
Re p = U s d p v~ l >^ 1, equation ([1]) becomes inaccurate because of nonlinearities that 
appear in the drag due to boundary layer separation. These are quantified with empirical 
correction factors, such as (3 = (1 + 0.15i?ep 687 ) [5] which appears in the drag model as: 



F D = 37Tfifd p U s = aU, 




F = aU s /3(U s ), 



(2) 
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Neither equation (1) nor equation (2) is valid for nonuniform flows. In the case of simple 
shear, some extensions are possible . Turbulence, however, is far more complicated, for 
the flow exerts time- varying stresses on the particle surface that are completely different 
from those due to boundary layer separation in uniform flow. The resultant force on a 
particle then includes a history term and an inertial term that depend on particle-fluid 
relative accelerations [61 [7J . 

Small particles often experience small relative acceleration, and see the turbulent 
stresses as uniform and linear fields. As a result equations ([I]) and (j2J) stay approximately 
true. The slip velocity in this case can be defined as the instantaneous difference 
between the particle's Lagrangian velocity u p l and the co-located Eulerian velocity of 
the fluid Uf, giving Use-l = u f ~ u pL- This definition works because the fluid velocity 
field can be considered 'unbroken' by the presence of the particles. Use-l is used in 
numerical models that perform Lagrangian tracking of infinitesimally small particles 
through Eulerian velocity fields [3]. The results have led to some significant advances 
in understanding, but also deserve a significant amount of caution, as discussed in [Hj. 

Larger particles experience nonuniform stress fields due to turbulence. There may 
be some large lengthscale at which the turbulent field can be filtered such that the 
particle responds to the filtered velocity field Uf according to the drag law: 

F d = aU s (3(U s ). (3) 

However, the particle will also experience drag due to the sub-filter-scale stresses. This 
additional drag will be relevant to processes such as particle collisions and fluid TKE 
modulation. The form of the sub-filter-scale terms is far from trivial, as expressions 
for history and inertial terms due to the unsteadiness are complex and their analytical 
expression has only been derived for the special cases of tiny particles moving with small 
slip velocities [9]. The work of Bagchi & Balachandar (TO] provides a good foundation 
for understanding the effect of particle size and turbulence intensity. In this work, the 
authors explored the case in which a particle is unmoved by any turbulent fluctuations, 
but rather stays fixed as a turbulent flow advects past at a constant speed. They 
use direct numerical simulation (DNS) to exactly determine drag on fixed particles 
in turbulent flow and compare the results to drag models. Particles have dimensions 
dp/ 7 ] — 1-53 — 9.59 and are subjected to meanflow with free-stream turbulence ranging 
from 10-25% (Re p =60-600). They find that drag models of the form of equation (j2J) work 
well for predicting the time-averaged forces for small particles. They also show that such 
drag models perform reasonably well for predicting instantaneous forces experienced by 
small particles, as long as the correct instantaneous slip velocity is used. 

As particle size and turbulence level is increased, Bagchi & Balachandar find 
that the neither drag models given by equations (j2J) or (j3J) are able to capture drag 
fluctuations observed in the DNS. Adding history terms derived by [7J does not improve 
the predictions, but rather introduces unphysical oscillations in the predicted drag. 
An obvious hurdle they face in testing the drag model is choosing a definition of 
instantaneous slip velocity. This is because the large particle size and the multiscale 
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nature of turbulence prevent an objective definition of local and undisturbed flow 
velocities. In this case drag coefficients become dependent on the choices of the volume 
around which averages are performed. Another problem is that the kernels used to 
compute history terms were derived assuming small Reynolds numbers and slowly 
varying slip- velocities, which is not valid when turbulence intensity is high and particles 
are large. 

These results indicate that we must reconsider both drag models and the definition 
of characteristic slip velocities for the case of large particles in turbulent flows. A 
possible approach is to model turbulent drag as a combination of steady and fluctuating 
terms. Earlier we defined U s as the difference between 'appropriate' fluid and particle 
velocities. To perform a Reynolds decomposition, we work with expectation values 
(■), which practically are approximated as averages over a sufficiently large statistically 
homogeneous region in space and/or a statistically steady period in time. We perform 
the Reynolds decomposition on instantaneous fluid velocities Uf and instantaneous 
particle center of mass velocities u p . The residuals u'* and u' p can be found from the 
instantaneous values it, — but we note that co- located values u' p — u'f are only 
clearly defined if (•) is replaced with a spatial filtering operation. For the remainder of 
the work presented here, (•) is sufficient for our purposes. Thus, performing a Reynolds 
decomposition on the slip velocity gives U s = ((w/) — {u p )) + (u'f — u' p ) = (u s ) + (u'f — u p ). 
We can use the two slip velocities to decompose drag into steady and fluctuating terms: 

F D = F + G, (4) 

F D = a(u s }f3((u s }) + G, (5) 

where the effect of the fluctuating component of the slip velocity {u't— u' p ) is included 
in G. The steady term has been shown to work relatively well in some cases [TUl [TT] . 
thus G can be considered either a correction term or a noise term [T2l [13]. It can be 
formulated in analogy to the standard drag law: 

G = a(u' f - u' p )P{u' f - u' p ) (6) 

Or given its own formulation 

G = f(u' f , u' p , d p , p f , [if, . . .) (7) 

G can be deterministic, but that carries the same problems as specifying Fjj. 
Alternatively, G can be considered a stochastic term. Because the idea of slip velocity 
is compelling to our intuition, and related to the idea of a relaxation time, we would 
like to identify a single velocity scale to use in / that will describe the stochastic "slip" 
between particle and fluid velocity fluctuations. 

3. Stochastic slip velocity 

Two obvious choices for the stochastic slip velocity V are: 

V^««?>-<«?» 1/a (8) 
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and 

V B = {(u' f -u' p ) 2 ) 1/2 . (9) 

We consider both of these definitions. One immediate advantage of Va is that 
it can be determined by measuring the fluid and particle phase velocity statistics 
independently. In contrast, Vb includes a covariance (uffu p f). An advantage of Vb is that 
it arises naturally from equation El being the variance of the instantaneous slip u'f — u' p . 
A third option, Vq = (u' 2 )' 1 ^ 2 — (u 12 )' 1 ^ 2 , is an algebraic variation on Va- We prefer Va 
over Vc because Va can be easily visualized as the difference between the power spectrum 
of particle velocity and the power spectrum of the fluid phase velocity. In addition, Va 
can be calculated equally well for laboratory, field, and numerical experiments. 

To understand the behavior of Va and Vb, we consider two extreme cases. A passive 
tracer, by definition, will have (u' 2 ) = (u 12 ) = (u'fU p ). In this case, both stochastic 
fluctuating slip velocities will be zero: Va = Vb = 0. The other extreme is a particle 
that does not respond to any of the fluctuating fluid velocities. This may be because 
the particle is anchored at a location, or because it moves "ballistically" through the 
flow with large inertia. By this definition, this extreme case has u' = 0, which makes 
{u'fU' p ) = as well. Thus in this extreme case both our definitions of slip velocity equal 
the magnitude of velocity fluctuations in the fluid phase: Va = Vb — (u' 2 ) 1 ^ 2 . Both of 
these limiting behaviors (totally passive and totally ballistic) yield stochastic fluctuating 
slip velocities that make good sense: passive particles do not slip relative to the fluid 
phase, while 100% of the turbulent fluctuations slip past ballistic particles. 

We can also use these limiting behaviors to provide a constraint that can help us 
prescribe models for {u'fU' p ). We start with the assumption that {u'fU p ) = j\{u' 2 ), (u' 2 )). 

Nondimensionalizing by (u' f 2 ) gives ," p s - = F( pf^r )- The limiting behaviors discussed 
above imply that F(0) = and F(l) = 1. The simplest function satisfying these 
constraints is linear: F(x) = x. This choice implies that the covariance equals the 
particle velocity variance, which makes the definitions for Va and Vb identical. Until 
proper data is available to support a more complex model for the covariance, we will 
assume that this simple model is acceptible. This leaves us with a single definition of 
stochastic slip velocity: V = Va ~ Vb- We analyze V using laboratory data below. 

4. Experiment 

Experiments are performed in a turbulent water tank, which uses two symmetric arrays 
of randomly firing synthetic jets to generate homogeneous isotropic turbulence with 
almost zero mean flow [14] . All measurements are taken within the region of isotropy and 
homogeneity, which is a cube roughly 30 cm on each side. In this region, the longitudinal 
integral length scale T x = 7.5 x 10~ 4 m, turbulent kinetic energy q 2 = 2.32 x 10~ 4 
m 2 s -2 , dissipation rate is e = 1.52 x 10~ 5 m 2 s~ 3 and the Taylor-scale Reynolds number 
R\ = 160. The facility is described in detail in [151 US]- 
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Table 1. Relevant particle parameters relative to turbulent scales. 



A suspension of large near-neutrally buoyant particles is studied at a volume 
fraction of 0.14%. Particles' density p p = 1020kg/m 3 at 20°C and they are manufactured 
by hand using injection molding. Particles are made of a transparent hydrogel solution, 
closely matched to the refractive index of water. Each particle contains hundreds of 
small (~ 40 /im) optical tracers, allowing velocities to be measured within each particle 
by tracking the tracers. 

Two particle shapes are considered here: spheres of diameter d s — 8 mm, and 
prolate ellipsoids whit polar and equatorial axes l e = 16 and d e — 8 mm, respectively. 
These particle sizes are within the inertial subrange of ambient turbulence, between the 
Kolomgorov microscale (77 = 0.38 mm) and the integral lengthscale (r^ = 75 mm), and 
of similar size to the Taylor microscale (A^ = 16 mm). More details are given in table 

m 

Although Stokes-based particle relaxation time is not useful for predicting the 
dynamics of particles of this size, we calculate it here as a point of reference. The 
Stokes-flow relaxation time of the spherical particles is Tpi = 3.64 s. The Stokes- 
fiow relaxation time of the ellipsoidal particles is computed using the expression for the 
average drag coefficient of randomly oriented prolate ellipsoids found in Clift et al. [5]: 



/ 2 T 

C el = 3vr4 Va ~ (10) 

log (a + va 2 — 1) 

where a = l e /d e is the particle aspect ratio. Consequently, the particle response time 
for our ellipsoids of a = 2 and d e = d s becomes: 

rW=rWa!55&±^S = i.B^) (11) 

where 7^ is based on the smaller (equatorial) diameter of the ellipsoid. Our ellipsoidal 
particles thus have Tpf = 5.46s. According to [5], we can define an equivalent diameter 
d p eq ^ for ellipsoidal particles as the diameter of a sphere having the same average drag 
coefficient as the ellipsoid. This can be determined from equation [10] to be d p eq ^ = 10.5 
mm. Equation ( fill offers an interesting insight on the effect of particle shape: particle 
response time is the ratio between inertial and drag force; while this ratio scales with 
particle size for spheres, this is not true for elongated particle. For example, increasing 
particle size by increasing the aspect ratio might actually decrease the response time 
and the equivalent diameter. 

Fluid phase velocity measurements are performed using standard 2D particle image 
velocimetry (PIV). These measurements use a single camera oriented perpendicular to 
a measurement volume of dimensions 4 cm x 3.4 cm x 0.1 cm centered in the tank. 
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Table 2. Second and fourth order statistics of particle and fluid phase velocity. 

The final size of the interrogation area is 32x32 pixels, corresponding to 0.8x0.8 mm in 
physical space. This size is on the order of 2r), thus spatial resolution effects are expected 
to be negligible. Statistics of the fluid phase are based on 510 independent velocity fields 
(corresponding to > 300 integral time scales) containing 102x86 datapoints each. 

We measure velocities within particles using stereoscopic PIV. The optical setup 
is two cameras (Imager PRO-X, 1600 x 1200 pixels, both fitted with a 50 mm Nikkor 
lens and Scheimpflug/tilt adapter) viewing a 1mm thick laser light sheet, capturing 
a measurement volume of 14 cm x 8 cm x 0.1 cm centered in the tank. The two 
cameras view the measurement area from opposite sides, each at an angle of 35 degrees 
relative to the laser's forward-scatter direction. From the velocities within each particle, 
we compute linear and angular velocities of the particles (see Appendix). Statistics 
of angular velocity are not presented here, but they are important when measuring 
large particles. For example, we have used them to control error in the linear velocity 
measurements (see Appendix). Statistics of particle velocity are obtained from 927 
and 1076 trajectories for spheres and ellipsoids, respectively. More details on the 
experimental setup and measurement techniques are given in |17j . 

5. Result 

The probability density functions (pdfs) of particles' linear velocities (u p ) are shown in 
figure [Th. For both spheres and ellipsoids, the pdfs are symmetric around u p = 0, as 
expected for homogeneous isotropic turbulence with zero mean flow. Table [2] shows that 
both particle types have velocity variance that is significantly smaller than that of the 
fluid phase. This is an immediate indication that particles do not behave as passive 
tracers. The variance of ellipsoidal particle velocities is significantly lower than that for 
spheres, implying that they follow the flow less faithfully than the spheres, likely due 
to their shape or to their longer relaxation timescale (ellipsoids' Stokes-based response 
time is 50% longer than spheres, as derived in section H]). 

Figures [T)d and [Ik compare the particle velocity pdfs to the fluid phase velocity. We 
include pdfs of the velocity field filtered at lengthscales corresponding to the spherical 
particle diameter ((w/) d ) and the ellipsoids' major axis ((u/);). Comparison with these 
demonstrates the nonlinear dynamics of particle-turbulent interactions, as discussed 
below. When particles show preferential concentration, they selectively sample the fluid 
velocity field [181 CHI EQ] . Based on the results of [21], we do not expect that our particles 
will behave in such a way, and thus we assume that they uniformly sample the velocity 
field. If not, the particle pdfs should be compared to fluid velocity statistics conditioned 
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Figure 1. a) Histograms of spherical (solid lines) and ellipsoidal (dashed lines) 
particles' velocity components, b) Comparison with fluid velocity pdfs. c) Normalized 
pdf showing hyper-gaussianity of particle velocity. The figure shows that particle 
velocity is not trivially related to velocity at any scale. 



on vorticity [22] . 

The shape of the normalized pdfs seen in figure [Tfc is of particular interest, because 
it reveals whether particles simply act as a linear filter of the velocity field, or if the 
interaction is more complex. From figure Hb we can see that while fluid velocity follows 
a nearly gaussian distribution, particles show a super-gaussian behavior, which is even 
more pronounced for ellipsoidal particles. We quantify this effect with fourth order 
moments of the pdfs, see table [2]). This indicates that the filtering relationship between 
particle and fluid velocity is nonlinear, as is discussed further in section [6j 

Applying the definition of slip velocity given in equation ([9]) we obtain that slip 
velocity for spherical particles is = 4 mm/s and for ellipsoidal particles is V^ e > = 6.8 
mm/s. The definition of stochastic slip velocity V can also be modified to so that it 
uses the variance of a spatially filtered fluid velocity field. The advantage of doing this 
is that it could explicitly represent the manner in which particles spatially filter the 
fluid velocity field. Doing so also provides an unambiguous definition for the covariance 
term in Vb- However, it is not yet known what filter best approximates the particles 
experience of the fluid velocity field. Indeed, if this were known, it would be an important 
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Table 3. Slip velocities and particle Reynolds numbers. 



step forward in modeling fluctuating drag on large particles in turbulence. Here, we use 
a simple boxcar filter, and find that the slip velocity does not change significantly with 
this filtering (see table [3]). This is likely because the variance is a large-scale statistic and 
particles are at the Taylor microscale. Filtering the fluid at progressively larger scales 
will eventually cause the slip velocity to become negative. Apparently by coincidence, 
we have found that this critical filter size for spherical particles in our flow is twice the 
sphere diameter. 

We can use the stochastic slip velocity V to define a Reynolds number Rey = 
Vd v v~ x . For the particles measured here, Rey is 32 for spheres and 69 for ellipsoids. We 
consider this result the first indication that Rey may indicate the degree of nonlinearity 
in the particle fluid interactions. This hypothesis is driven by two observations. First, as 
Rey approaches zero, the fluid particles follow the flow as passive tracers, and thus their 
motion is a linear function of the fluid phase velocity. Second, the particles measured 
here have Rey > 1, demonstrate a nonlinear response to the fluid velocity field (as 
indicated by the pdfs), and the strength of this nonlinearity is correlated with the 
Reynolds number (i.e. ellipsoids have a larger Rey and exhibit greater nonlinearity). 
Of course, this correlation is based on only two points in parameter space, and our 
assumptions about how Rey should indicate drag linearity are based primarily on 
analogy with the particle Reynolds number in Stokes flow. Thus more analysis is 
needed to understand the meaning of Rey, which we expect will be closely tied to 
the development of stochastic drag models employing V. 

6. Discussion 

The kurtosis and scaled pdfs of particle velocities show that large particles like the 
ones used here exhibit super-gaussian behavior. The significance of super-gaussian 
behavior is that it indicates intermittency. Intermittent behavior is typically seen 
in small-scale turbulent processes [231 I2H 125] and thus the intermittency in particle 
dynamics may help identify links between particle motion and fluid flow. Super- 
gaussianity of large spherical particles was also reported by Qureshi et al. [I]. In their 
experiments, intermittent behavior persists for particle sizes up to the largest turbulent 
scales, while fluid velocity increments at the scale of large particles are nearly gaussian. 
This indicates a non-trivial link between particle and fluid velocity. They suggest that 
particle acceleration statistics are related to statistics of the eulerian pressure increment, 
which indeed is a complex function of the velocity field. 

Particle velocity intermittency also appears to depend on particle shape. Our 
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data show that intermittent behavior is more pronounced for ellipsoidal particles than 
for spheres. We attribute this difference at least partially to shape, because the 
ellipsoids considered here have an equivalent diameter similar to the spheres' diameter 
(d^ = 1.3d p ^). Unlike spheres, ellipsoids have orientation-dependent drag. That is, 
in equation 1, the drag must also include orientation <fi: = f(Re p ,4>(t),...). As 
orientation fluctuates, so will the drag coefficient. 

The difference between the scaled particle and fluid pdfs in figure [lb might help 
point the way to an improved use of Stokes number for large particles in turbulent 
flow. The Stokes number is successful in predicting particle dynamics and turbulence 
modulation for small particles in turbulent flow [31 El EE]. When particles become 
large, then the Stokes number fails in two ways. First, the assumptions inherent in its 
derivation become invalid. Second, it fails to collapse dynamics across a wide range of 
flows. Lucci et al. (27] examined this failure, showing that large particles with identical 
values of Stokes number had very different behaviors. Similar problems with Stokes 
number are also seen for dense particles, but due to different physical effects, and are 
discussed in [28] . 

The definition of Stokes number uses an exponential particle relaxation timescale, 
which follows from a drag law that is linear in U s . If the steady drag F is similarly linear 
in (u s ), then an analogous result will follow. This is essentially what Xu & Bodenschatz 
[2Tj propose. They formulate a Stokes number that is the ratio of to a filtered turbulent 
timescale and particles' response time to a filtered velocity field. 

The fact that particle velocities do not follow the filtered velocity field (as seen 
in fig. [1]) implies that either F is nonlinear, G is non- negligible, or both. As a result 
we can predict that the filtered Stokes number will not effectively collapse dynamics 
that depend on the sub-filter scale motions. Notable effects that will depend on these 
motions are particle clustering and turbulence modulation. 

At this stage, we see three main uses for our definition of stochastic slip velocity 

1 /2 

V = \{u'j) — (Up)) . The first is to quantify the importance of fluctuating drag G. 
When V — 0, the particle acts as a tracer and thus G — 0. When V is nonzero, we 
can compare it to the steady slip velocity (u s ), to understand the relative importance 
of steady and fluctuating slip. As these slips enter into the steady and fluctuating drag 
(F and G) respectively, we can infer that FjG is a function of (u s )/V. Rather than 
comparing the slip velocities directly, we can also compare their equivalent Reynolds 
numbers: Rey = ((uf) — (u p )) 1 ^d p v~ 1 and Re s = ((uf) — (u p )) d v v~ l . 

The second major use for V is as an input for formulating G. There are many 
possible expressions for G, but we expect they all will include a comparison of (u"j) and 
(Up). We suggest that the difference V is preferable to the ratio (uj)^^) -1 because its 
form is analogous to the slip velocity that shows up in the standard and steady drag 
laws. Future work, primarily DNS, will be needed before the exact relationship between 
G and V can be determined. 

The third major use for V is in predicting some of the key features of particle- 
laden turbulent suspensions. There is currently no method for predicting turbulence 
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modulation (such as TKE amplification or attenuation) and particle clustering (i.e. 
preferential concentration) for suspensions of large particles in turbulence. For small 
particles, the Stokes number has been used to predict these effects, but it fails for large 
particles as discussed above. We expect that whatever non-dimensional parameters 
predict clustering and TKE modulation for large particles will also control the stochastic 
slip velocity V. If this is true, then V could serve as a midway predictor for clustering 
and TKE modulation. This would serve two purposes: V may be easier to measure 
than the clustering and TKE modulation, and the relationship between G and V, once 
found, would guide the search for a predictor of clustering and TKE modulation. 

7. Conclusions 

The standard drag model fails to accurately describe forces on large particles in turbulent 
flows. We can replace the standard drag with the sum of steady and fluctuating drag 
terms, but the fluctuating drag term has not yet been parameterized. Herein, we discuss 
the definition of a stochastic slip velocity. This can serve as an input to fluctuating drag 
models, and also serve other key analyses in turbulent particle suspension dynamics. 
Specifically, it allows us to determine a particle Reynolds number for neutrally buoyant 
particles, and contributes to the extension of the Stokes number. This is essential for 
efforts to accurately predict the dynamics of large particles in turbulent flows. 

We discuss two options for the stochastic slip velocity, and demonstrate their 
approximate equivalence. We then measure the stochastic slip velocity using laboratory 
data collected on large ellipsoidal and spherical particles in homogeneous isotropic 
turbulence. The results suggest that particle Reynolds number based on the stochastic 
slip is an effective predictor of nonlinear fluid-particle interactions. Specifically, both 
particles show velocity distributions with greater kurtosis than the fluid phase velocity 
distribution. Comparing the ellipsoids and spheres, we see that the increase in kurtosis 
is correlated with an increase in particle Reynolds number. 

A key next step is to extend this argument to filtered, rather than expectation, 
velocities. This will allow us to quantify the covariance between velocity and fluid 
fluctuations, and thus evaluate the difference between our two proposed definitions of 
stochastic slip velocity (Va and Vb)- In addition, we can evaluate rotational slip by 
using within-particle velocities to extract particle rotation measurements. 

Appendix 

We use stereoscopic PIV to measure the velocities Um and at two points within 
the particle: Xm and Xjy- From these measurements the angular velocity O can be 
determined from the equation of solid body rotation: 



Since PIV measurements are in a single x — y plane, equation [12] is over-determined in 
Q z and under-determined in Q x and Q y . By including a third point P in addition to M 




(12) 
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and N, all three components of f2 can be determined. 

Because PIV provides many more than 3 velocities within each particle, we can 
improve the precision and accuracy of ft measurements. We do this by calculating a 
value of O using each possible set of three points M, N and P, and considering the 
distribution of Q estimates. 

The location of the particle center of gravity is not explicitly measured by PIV, 
but the velocity there (V cg ) can be determined as follows. The equation of solid body 
rotation gives: 

V cg = V cf - n x r, (13) 

where V c f is the velocity at the center of the particle cross-section formed by the 
intersection of the particle and the PIV measurement plane, fl is the angular velocity 
as calculated above, and r = (r x ,r y ,r z ) is the distance between the center of gravity of 
the particle and the center of the particle cross-section. Note that for spherical particles 
r x = r y — 0. From equation (fTSj) . we can express the variance of V cg as: 

(V cg 2 ) = (V cf 2 ) - <(n x if) - (2V cg (fi x r)>. (14) 

It is reasonable to assume that the last term on the RHS of equation [T41 vanishes because 
the velocity of the center of gravity and the velocity induced by particle rotation at the 
measurement location are uncorrelated. We therefore are left with the expression: 

(v cg 2 ) = (v cf 2 ) - ((n x if). (is) 

This indicates that the variance of the velocity at the particle's center of gravity is the 
variance of the velocity measured at the center of a randomly oriented plane intersecting 
the particle minus the 'rotation-induced noise term ((fi x r) 2 ) that is due to the velocity 
induced by particle rotation at the measurement point. 

To quantify the importance of the rotation-induced noise term, we expand V r i n = 
f2 x r in its components. Because spherical particles have r x = r y = 0, the noise 
term becomes (u r i n ,v rin ,w r i n ) = (Q y r z , —Q x r z ,0). Therefore, we can immediately 
conclude that the variance of w c f equals the variance of w cg . The variance of the other 
two components can be expanded by noting that the vector f2 is independent of r z . 
Therefore we can express the variance of the products as the product of the variances: 
((Qir z ) 2 ) = (Q 2 )(r 2 ). The variance of the distance r z of a point chosen randomly within 
the volume of a spherical particle of radius i? is (r 2 ) = |i? 2 . Therefore (u 2 rin ) = |_R 2 (f2 2 ), 
and (v 2 J = §i? 2 (fi 2 >. 

In our measurements, the diameter of spherical particles is R = 4 x 10 -3 m, and 
the angular velocity measurements show that (Q 2 ) = O(10 _1 ). Therefore the noise term 
is O(10 -7 ). Comparing this to (u 2 ^), which is O(10~ 4 ), we conclude that the rotation- 
induced noise is negligible in our measurements of spherical particles velocity variance. 
For ellipsoidal particles, the expression for V r i n becomes much more complicated because 
r x and r y are non-zero and the expression for the radius variance is different than that 
for spheres. However, our measurements of O for ellipsoids are very similar to the 
results for spheres. Given this, and the fact that the ellipsoids are very similar size to 
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the spheres, we expect the order of magnitude of the noise term to be similar for spheres 
and ellipsoids. As a result, we neglect the rotation-induced noise when calculating the 
variance of ellipsoid velocities. Should these terms be non-negligible in other cases, they 
can be quantified using the method presented here. 
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